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We apply the large-deviation method to study trajectories in dissipative quantum systems. We 
show that in the long time limit the statistics of quantum jumps can be understood from thermo- 
dynamic arguments by exploiting the analogy between large-deviation and free-energy functions. 
This approach is particularly useful for uncovering properties of rare dissipative trajectories. We 
also prove, via an explicit quantum mapping, that rare trajectories of one system can be realized as 
typical trajectories of an alternative system. 



Introduction. Equilibrium statistical mechanics provides 
the tools to study equilibrium phases and phase changes 
in many body systems pp. Thermodynamic phases are 
characterized by average values of thermodynamic ob- 
servables, such as volume in a liquid or magnetization in 
a magnet, which are controlled by conjugate fields, such 
as pressure or magnetic field. Non-analyticitics in free- 
energies correspond to phase transition points, and the 
proximity to a phase transition manifests in large and 
rare fluctuations of observables around their thermody- 
namic values pp. 

An analogous perspective can be adopted for the study 
of dynamical phases in non-equilibrium systems by ap- 
plying the large-deviation (LD) method [2J. The LD for- 
malism allows to treat ensembles of trajectories, classified 
by dynamical order parameters or their conjugate fields, 
in the same way that equilibrium statistical mechanics 
treats ensembles of configurations. Important properties 
of classical non-equilibrium systems can be uncovered by 
exploiting this analogy [2, 3, 4 , such as the existence of 
"space-time" phase transitions in glassy systems [5]. 

In this Letter we apply the LD method to quantum 
non-equilibrium systems. This approach reveals impor- 
tant properties of ensembles of trajectories of quantum 
systems that undergo quantum jumps in some form, such 
as driven quantum systems weakly coupled to a thermal 
bath [51 \7\. We show that one can observe features of dy- 
namical crossovers and dynamical phase transitions even 
in quantum systems with only a few degrees of freedom, 
and illustrate our ideas with three simple examples: (i) a 
driven 2-lcvcl system, where the LD approach allows to 
identify a scale invariancc point in the ensemble of trajec- 
tories of emitted photons; (ii) a blinking 3-level system 
(or electron shelving problem), where we argue that in- 
tcrmittency in photon count is related to a crossover be- 
tween distinct dynamical phases; and (iii) a micromaser, 
where static bistability leads to a first-order phase tran- 
sition in the ensemble of trajectories. We also establish 
a mapping between two dynamical systems, where typ- 
ical trajectories of one are the rare trajectories of the 
other. This method is particularly useful for generating 
rare trajectories which otherwise are highly suppressed. 

Formalism. We consider a quantum system weakly cou- 
pled to a reservoir in the Markovian regime. The non- 



unitary evolution of its density matrix p(t) is described 
by a so-called Master equation jSJ E] , 
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where L M and (p, — 1, . . . , N-^) are the Lindblad op- 
erators [EJ |9], {•,•} stands for anticommutator, and we 
have set h = 1. We are interested in the time record of 
projection events due to one (or more) of the Lindblad 
operators, such as that of emitted photons, which we as- 
sume are detected with 100% efficiency. Such record is a 
particular quantum jump trajectory of the system [9]. 
The probability Pt(K) to observe K events after time t 
is given by Pt(K) = Tr [p^(t)l, where p( K \t) is a re- 
duced density matrix obtained by the projection of the 
full density matrix onto the subspace of K events, e.g., 
the subspace containing K photons [6]. For large times 
Pt(K) acquires a LD form: 
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The "large-deviation" function ip(k) (k = K/t) contains 
all information about the probability of K at long times. 
Alternatively, we can describe the statistics of K via the 
generating function, which also has a LD form, 



Zt(s) 



OO 

£ 

K=0 



Pt(K)e 



~sK 



J8(s 



(3) 



The LD functions f(k) and 9(s) are to trajectories what 
entropy density and free-energy density are to configu- 
rations in equilibrium statistical mechanics, with s be- 
ing the conjugate field to the dynamical order parame- 
ter K. The two are related by a Legendre transform, 
9{s) = — mint [<p(k) + ks], and the function 9(s) has the 
convexity properties of (minus) a free-energy. Moreover, 
anomalous dependence of 9(s) on s indicates non-trivial 
fluctuation properties of dynamical trajectories. In par- 
ticular, singularities in 9(s) correspond to dynamical (or 
space-time [5]) phase transitions. It is this anomalous 
and phase-transition behavior that we uncover below by 
calculating 9(s) for simple driven quantum systems. 

The p^ K ' (t) obey a set of equations 6. which is uncou- 
pled by the Laplace transform |3]). That is, the equation 
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FIG. 1: (A) Laser driven 2-level system coupled to a T = 
bath. (B) Large-deviation function 8(s) of number of emitted 
photons K. Dynamical trajectories go from more active to less 
active as s, the conjugate field to K, is increased, as shown by 
the average photon rate k(s) = {K) s /t — —8'(s). The Mandel 
parameter Q(s) — — 2/3 for all s, indicating that for k = 4f2 
trajectories display a form of scale invariance. (C) The photon 
count probability is obtained from Q by a Legendre trans- 
form: P t (K) w e - tviK/t) with tp{k) = 3[k\n(k/k )-(k-k )]. 
It is a v = 3 Conway-Maxwell-Poisson distribution [14] , 
Pt(K) oc [Poisson^; t)] 3 . (D) Representative trajectories 
from sub-ensembles with different average k. 



for Ps (t) ee EK= P {K) (t)e- sK reads f t p s {t) = W s ( Ps ), 
where the super-operator W s is 

W s (p) = -i[H,p] + e- s L lP L\ 
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L\ is the Lindblad operator which produces the quan- 
tum jumps we are counting in K . The equation dtp = 
W s (p) is sometimes called the generalized quantum mas- 
ter equation [TOl HI]- The operator W s is analogous 
to the Lebowitz-Spohn operator [T^] of classical non- 
equilibrium dynamics. Physical dynamics takes place 
at s = [here W s = W, see Eqs. ( [Tp] )]. For s ^ 0, 
<9(p = W s (p) describes a time evolution whose unfolding 
[5J HI [13] generates an ensemble of trajectories biased by 
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see Eq. (|3h . We call this the s-ensemble |5J. 



In analogy to the classical case we assume that the LD 
function 8(s) is given by the largest real eigenvalue of W s 
|12j . This assumption, corroborated below, reduces the 
problem of calculating the "partition sum' ' of Eq. ([3]) to 
an eigenvalue problem. 

(i) 2-level system. Consider a 2-level system, Fig. [TJ\, 
driven by a resonant laser in contact with a zero temper- 
ature bath [9j. When the observable K is the number of 
emitted photons the generalized master operator is 

W s (p) = —ifl[a + a\p] + e~ s n apa) — — (a) ap + pa* a) , 



where a and a* are the lowering and raising operators, 
|0)(1| and |1)(0|, respectively, ft is the Rabi frequency, 
and k is the decay rate. We consider the specific choice 
k = 4f2, which is interesting for reasons we discuss below. 
Here the LD function takes the simple form 
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which is shown in Fig.[T}3. R vanishes at s = 0. This is a 
statement of conservation of probability: Wo reduces to 
the master operator of Eq. Q which leaves Tr[p] invari- 
ant. Derivatives of 6(s) give moments of the photon num- 
ber distribution. In particular, the average number of 
emitted photons is ko = (K)/t = — 0'(0), and the Mandel 
parameter, Q = ((K 2 ) - (K) 2 )/(K) — 1 = -0"(Q)/0'(O). 
The LD function around s — encodes the information 
about fluctuations of typical trajectories |4j [10] . 

Away from s = 0, 8(s) encodes information about rare 
trajectories. Consider the s-dependent average photon 
number (per unit time), 
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This expression is the average of K/t where the probabil- 
ity of trajectories is biased by the factor e~ sK . Pursuing 
a thermodynamic analogy, think of K and s as volume 
and pressure. Increasing/dereasing pressure leads to a 
smaller/larger average specific volume, i.e. by control- 
ling pressure we obtain a denser or less dense system. 
Something analogous occurs here in the dynamics: s > 
corresponds to trajectories with fc(s) < fco, i.e. less ac- 
tive than typical, while s < corresponds to trajectories 
with k(s) > k , i.e. more active than typical, Fig. [lj3. 

We can also define an s-dependent Mandel parame- 
ter, Q(s) = {{K 2 ) s -{K)l)/{K) s ~l = -6"(s)/6>( S )- l, 
which measures the bunching or anti-bunching properties 
of trajectories with a fixed average photon number tk(s). 
For the specific case of n = 40 we have k(s) — 20e~ s / 3 /3 
(i.e. trajectories go from more to less active as s is in- 
creased from negative to positive), but Q(s) — —2/3 for 
alls, Fig.[l£. This result is surprising. We expect photon 
emissions to be anti-bunched [9] , but an s-indepcndent Q 
indicates that all sub-ensembles of trajectories, no matter 
how active or inactive, have the same fluctuation prop- 
erties of typical trajectories: trajectories would look the 
same if rescaled by their average emission rate. Hence 
k = 4D, is a "special point" in parameter space where the 
dynamics displays trajectory scale invariance. Note that 
this occurs while all correlation times remain finite. 
(ii) 3-level system. Consider now a 3-lcvcl system like 
the one of Fig. [2]A, driven by two resonant lasers on the 
|0)-|1) and |0)-|2) lines with Rabi frequencies fii and O2, 
respectively. Level |1) decays to |0) with rate K\. We are 
interested in the statistics of the number K of photons 
emitted. When fl\ 3> O2 typical photon emission trajec- 
tories are intermittent, displaying "bright" and "dark" 
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FIG. 2: (A) Laser driven 3-level system. Here Ki = 4f2i and 
0,2 = fii/10. (B,C) The LD function 8(s) and dynamical or- 
der parameter k(s) display crossover behavior near s = be- 
tween active and inactive dynamical regimes. The active side 
is anti-bunched, Q < 0. The inactive side is non-fluctuating 
Q = 0. The peak in Q near s — signals the dynamical 
crossover. (D) The fat tail for k < ko in Pt(K) is a mani- 
festation of the inactive regime; the thin tail for k > ko is a 
manifestation of the active regime. (E) Representative trajec- 
tories from inactive and active sub-ensembles. At s = there 
is (mesoscopic) coexistence of the two dynamical regimes and 
typical trajectories are intermittent or "blinking" . 



periods [7J[TS]. In this case quantum jumps can become 
evident on macroscopic timescales [To] . 

The generalized master operator W s is of the form Q , 
with H — 5Z J=1 ^j( a j + a ]), where a,j = \0)(j\ and dj = 
|j)(0|, and only one set of Lindblad terms, N-^ — 1, with 
L\ = y/~Kiai. The LD function 9(s) is obtained from W s 
by direct diagonalization. It is shown in Fig. [2j3 for n\ — 
4f2i and ft 2 = fii/10. The difference with the 2-level 
case is striking. For s < (i.e. trajectories more active 
than typical) 9(s) follow the LD function of the 2-level 
problem. Close to s = 0, however, 9(s) leaves the 2-level 
curve and approaches a constant, 9(s ^> 0) ~ — reif2§- 

Within our thermodynamic analogy this indicates a 
rapid crossover between two distinct dynamical phases 
as we cross s — 0. Figure [2p shows the corresponding 
change in k(s). The active side is that of s < 0, and 
trajectories have large K. The inactive side is s > 0, and 
trajectories have small K. The active phase is that of the 
2-level system |0), |1) where photon emission is plentiful. 
In the inactive phase the atom predominantly occupies 
the 1 2) state and photon emission is scarce. The crossover 
in k(s) is reminiscent of a (smoothed) dynamical first- 
order transition, such as that seen in the trajectories of 
certain glassy systems [5]. The dynamical crossover is 
also apparent in Mandel parameter, Fig. |2p. The active 
phase is antibunched, Q(s «0) = —2/3, while the inac- 
tive phase does not fluctuate, Q(s ^> 0) = 0. The peak in 
Q(s) around s — is a signature of the crossover between 
phases: here fluctuations are maximal as trajectories are 



(mesoscopic, i.e. finite time) mixtures of the two coexist- 
ing phases. Typical trajectories correspond to s = 0, but 
the crossover structure of the LD function 6(s) has an 
effect on the tails of the distribution P t (K) , as shown in 
Fig. [2j3. It has a fat tail for k < (K)/t [originating from 
9 < for s ^> 0], and a thin tail for large k [originating 
from 9 Ri 02-level for s <C 0]. 

(Hi) Micromaser. We now consider the problem of a mi- 
cromaser |17j . a resonant cavity coupled to a finite tem- 
perature bath and pumped by excited two level atoms 
which are sent into the cavity with a constant rate, Fig. 
|3}\. The atoms only interact with a single mode of the 
cavity. The cavity reaches a steady state which is sen- 
sitive on pump rate and atom-cavity coupling. In par- 
ticular, the steady state of the cavity can change from 
unimodal to bimodal [17] . We now show that this static 
bistability has an associated dynamic bistability. 

Our dynamical order parameter K is now the num- 
ber of atoms which leave the cavity and are in the 
ground state. The super-operator W s Q follows from 
the Lindblad master equation for the cavity after trac- 
ing out the atom and the thermal bath [T7J. There 
are four sets of Lindblad operators, 2Vl = 4, two from 

the atom-cavity interaction, L\ = Jr — v , — r — '-a and 

V flflT 

L2 = \/rcos ^4>V aotj , and two from the cavity-bath in- 
teraction, L 3 = y^a and L 4 = s/Xa) . Here a, are 
the raising/lowering operators of the cavity mode, r is 
the atom beam rate, k and A are the thermal relaxation 
and excitation rates, and <f> encodes the atom-cavity in- 
teraction |17] . Events are recorded when quantum jumps 
under the action of L± occur. 

The LD function 9(s) can be obtained by assuming 
that the corresponding eigenmatrix r s of W s (see below) 
is diagonal in a^a. It is shown in Figs. [3p,C for two 
values of the "pump parameter" a = <\>\Jr /(k — A) |17j . 
For a = 2ir the stationary state of the cavity is close to 
being bistable, undergoing a sudden change from a low 
average photon occupation (N) at a < 2ir to a large (N) 
at a > 2-7T 17J. In this case the LD function is singular at 
s = 0, and k(s) has a discontinuous jump, Fig. [3p. This 
is a dynamic, or space-time, phase transition 5 . It is 
first-order because the order parameter, k, changes dis- 
continuously. The phase transition is between an active 
phase at s < and an inactive (or less active) one at 
s > 0. Fig. [3p shows that the active phase corresponds 
to that of large (N) and the inactive one to small (N). 
The transition point is s = so that normal dynamics 
occurs under dynamic phase-coexistence. The dynamical 
transition remains even far from static bistability, but 
the transition point moves away from s = 0, Fig.[3}3. 

Mapping of rare trajectories to typical ones. The LD 
function 9(s) encodes properties of rare quantum trajec- 
tories, and by exploiting the analogy with thermodynam- 
ics we can describe sub-ensembles of trajectories as dy- 
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FIG. 3: Dynamical phase transition in the micromaser. (A) 
Cavity mode driven by pumped atoms and interacting with 
thermal bath. (B,C) LD function 8(s) for the number of 
atomic transitions, K. When the cavity is close to static bista- 
bility, a = 2-7T, the LD function has a first-order singularity at 
s — 0. There are two distinct dynamical phases, a more active 
one with large K, and a less active one with small K. Typi- 
cal trajectories are at coexistence between these phases. The 
dynamical transition is still present far from static bistability, 
a — 1.2n, but the transition point is at s < 0, i.e. dynami- 
cal coexistence will be only manifest in rare trajectories. (D) 
Cavity photon distribution in active and inactive phases, and 
at coexistence (i.e. stationary density matrix). 



The explicit construction of the trace-conserving sys- 
tem gives interesting insights into the structure of tra- 
jectories away from s = [18]: (i) For the 2-level system 
above we have H = e~ s ^ 3 fl(a + aJ) and L\ — e~ s ^y/~K,a. 
This is of the same form as the s — problem with 
all rates multiplied by e~ s / 3 . That is, rare trajectories 
(s ^ 0) are typical trajectories of the same system but 
with time rescaled as t — > e s ^ 3 t. This confirms k = 40 
as a special symmetry point of the 2-level problem, (ii) 
For the 3-level system above on the inactive side, s > 0, 
the mapped system is another 3-level problem with an 
additional strong laser on the |1)-|2) line. This coupling 
introduces an effective detuning for the laser on the |0)- 
|1) transition, suppressing the excitation of state |1) and 
subsequent photon emission [18] . 

The "statistical mechanics of trajectories" method pre- 
sented here seems to be a useful framework to study non- 
equilibrium dynamics in open quantum systems. Already 
for very simple problems it reveals unanticipated rich- 
ness. We expect it to be even more fruitful in the study 
of non-equilibrium quantum many-body systems. 

We thank Andrew Armour and Markus Miiller for 
helpful discussions. 



namical or space-time phases [S]. d t p = W s (p) however is 
not a physical time evolution, but we can show |18j that 
there is an alternative trace-preserving evolution which 
generates the same s-ensemble. Thus, rare trajectories 
in one system correspond to typical trajectories of a re- 
lated system, and crossovers or transitions controlled by 
s, such as the ones discussed above, can be realized as 
transitions controlled by physical parameters. 

The super-operator W s has the LD function 6(s) as 
its largest real eigenvalue, with "right" and "left" Her- 
mitian eigenmatrices r s and l g , respectively |19j . These 
cigenmatrices obey W(r s ) = 8(s)r s and (l s )W s = l s 9(s) 
and we normalize them such that Tr[r s ] = 1 and Tr[i s ] = 
Tr[I]. Given a matrix p s (t) which evolves according to 
dtp s = W s ( j o s ), there is an associated density matrix 
p(t) = l X J 2 p s (t)ll^ 2 /Tr[l s p s (t)], whose corresponding evo- 
lution, dtp = VV(/5), is of the Lindblad form (jlj, with the 
following Hamiltonian and Lindblad operators: 

H = \l- 1/2 ({H,l s }+ % -[LlL^l s ]\ (6) 



L, = 



<Ve~ s/2 + (1 - 



l\' 2 L^I\ (7) 



This dynamics is trace-preserving, (I)KV = 0, and the 
set of trajectories of quantum jumps due to L\ coincides 
with the s-ensemble of W s . The tilde process is that of 
a physical dynamics. Its typical trajectories correspond 
to rare trajectories of the original process W [20 . 
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